AEM Accepted Manuscript Posted Online 14 July 2017 

Appl. Environ. Microbiol. doi:10.1128/AEM.01326-17 

Copyright © 2017 American Society for Microbiology. All Rights Reserved. 


1 Title: 

2 Novel Alphacoronaviruses and Paramyxoviruses co-circulate with Type 1 and 

3 SARS-related Betacoronaviruses in synanthropic bats in Luxembourg. 

4 

5 

6 Maude Pauly a# , Jacques B. Pir b , Catherine Loesch 3 , Aurelie Sausy 3 , Chantal J. 

7 Snoeck 3 , Judith M. Hubschen 3 , Claude P. Muller 3 

8 

9 Infectious Diseases Research Unit, Department of Infection and Immunity, 

10 Luxembourg Institute of Health, Esch-sur-Alzette, Luxembourg 3 , Musee national 
n d’histoire naturelle, Centre de recherche scientifique - section Zoologie des 

12 vertebres, Luxembourg, Luxembourg 13 

13 

14 Running Title: Viruses circulating among bats in Luxembourg 

is #Address correspondence to Maude Pauly, maude.pauly@lih.lu 
16 M.P. and J.B.P. contributed equally to this work 


l 


Downloaded from http://aem.asm.org/ on July 18, 2017 by UNIV OF WINDSOR 




<D 

C 


c 

o 


“O 

0 


to 

O 

CL. 


U 

to 

3 

C 

O 


|D 

0 


CL 

0 

U 

u 

< 


o &5 

1-8 

LU _q 

“O o 


“O 

<D 


Q_ 

Q_ 

< 


Q_ 

Q_ 

< 


17 ABSTRACT: 248 words 

is Several infectious disease outbreaks with high mortality in humans have been 

19 attributed to viruses that are thought to have evolved from bat viruses. In this study 

20 from Luxembourg, the genetic diversity and epidemiology of paramyxoviruses and 

21 coronaviruses shed by the bat species Rhinolophus ferrumequinum and Myotis 

22 emarginatus was evaluated. Faeces collection (n= 624) was performed longitudinally 

23 in a mixed-species colony in 2015 and 2016. In addition, faeces (n= 254) were 

24 collected cross-sectionally from six Myotis emarginatus colonies in 2016. Using 

25 degenerate primers in a nested format, an overall prevalence of 1.1% (10/878) and 

26 4.9% (43/878) was determined for paramyxoviruses and coronaviruses. Sequences 

27 of the partial RNA-dependent RNA polymerase and spike glycoprotein genes of 

28 coronaviruses, as well as of the partial L-gene of paramyxoviruses were obtained. 

29 Novel paramyxovirus and Alphacoronavirus strains were identified in different Myotis 

30 emarginatus colonies, and SARS-related Betacoronavirus strains were shed by 

31 Rhinolophus ferrumequinum. Logistic regression revealed that shedding of 

32 Alphacoronavirus was highest in July (OR: 2.8, p<0.01), probably due to 

33 periparturient stress. Phylogenetic analyses point at close virus-host coevolution and 

34 the high genetic similarity of the study strains suggest that the Myotis emarginatus 

35 colonies in Luxembourg are socially connected. Most interestingly, we show that 

36 bats also host type 1 Betacoronavirus strains. The high similarity of the spike gene 

37 sequences of these viruses with mammalian type 1 Betacoronavirus strains may be 

38 of concern. Both the SARS-related and type 1 Betacoronavirus strains detected in 

39 bats in Luxembourg may cross the species barrier after a host adaptation process. 
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41 IMPORTANCE: 144 words 

42 Bats are a natural reservoir of a number of zoonotic pathogens. Several severe 

43 outbreaks in humans (e.g. Nipah virus outbreak in Malaysia in 1998, and the almost 

44 global spread of Severe Acute Respiratory Syndrome in 2003) were caused by bat- 

45 borne viruses that were transmitted to humans mostly after virus adaptation (e.g. in 

46 intermediate animal hosts). Despite indigenousness of bat species that host viruses 

47 with suspected zoonotic potential and despite zoonotic transmission of European Bat 

48 Lyssavirus type 1 in Luxembourg, knowledge about the diversity and epidemiology of 

49 bat viruses remains limited in this country. Moreover, in contrast to other European 

50 countries, bat viruses are currently not included in the national surveillance activities 

51 of this land-locked country. We suggest that this gap in disease surveillance should 

52 be addressed as we show here that synanthropic bats host viruses that may be able 

53 to overcome the species barrier. 

54 INTRODUCTION 

55 Their ability to fly long distances and their longevity enable bats (Chiroptera) to 

56 spread viruses across time and space. Large colony sizes, close social interactions 

57 and co-roosting of different bat species favour intra- and interspecies transmission of 

58 viruses (1). Moreover and above all, low pathogenicity of viruses and viral 

59 persistence in bats are indicative of ancient co-speciation between bats and different 

60 virus families [e.g. Paramyxoviridae and Coronaviridae (2-5)]. It has been suggested 

61 that most human coronaviruses (CoV) evolved from bat counterparts (5-7). For 

62 instance Severe Acute Respiratory Syndrome (SARS) (8) and Middle East 

63 Respiratory Syndrome (MERS) CoV (9, 10), but also Nipah and Hendra 

64 paramyxoviruses (PV) (11, 12) originated in bats and caused severe outbreaks in 

65 humans. While for some viruses viral adaptation processes in intermediate animal 
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66 hosts were presumably required before zoonotic transmission (e.g. (9, 13)), direct 

67 transmission of Nipah virus between bats and humans repeatedly occurred in 

68 Bangladesh (14). The spike glycoprotein of several bat CoV strains share features 

69 with human strains which were critical for bat-to-human transmission events (15). In 

70 particular, the receptor-binding domain of the spike gene determines the host range 

71 and tissue tropism of CoV (16-18). Nevertheless, the risk of a zoonotic infection with 

72 bat viruses is low in humans as direct contacts with bat excretions are rare (19, 20). 

73 In addition, the risk can be monitored by virus surveillance in synanthropic bats (20), 

74 such as vespertilionid (e.g. M. emarginatus) and rhinolophid bats (e.g. R. 

75 ferrumequinum) that have been shown to host a number of viruses with zoonotic 

76 potential (5, 21-23). 

77 In Western and Central Europe, M. emarginatus and R. ferrumequinum are 

78 endangered (48, 49) due to the on-going habitat fragmentation (24). After hibernating 

79 in underground sites, female R. ferrumequinum return in March to their natal 

80 colonies, while female M. emarginatus follow only in May (25-27). They form 

81 matrilineal maternity colonies in attics and barns (25-27). Around mid-June, females 

82 give birth to a single pup. Intra-lineage polygyny is common for R. ferrumequinum 

83 (28, 29) and extra-colony mating of R. ferrumequinum and M. emarginatus occurs 

84 during the swarming of the males, between September and October (30, 31). 

85 Despite a growing interest in these animals as hosts of emerging viruses, the 

86 knowledge about bat viruses in Luxembourg remains limited. In a single study 

87 European Bat Lyssavirus type 1 was isolated and the risk of zoonotic transmission in 

88 the country shown (32). 
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89 Here, we report the shedding of PV and CoV by R. ferrumequinum and M. 

90 emarginatus, two sympatric and synanthropic bat species. Virus diversity and 

91 prevalence was assessed in six nursing colonies of M. emarginatus in a cross- 

92 sectional manner. In addition, we investigated seasonal patterns of both viruses in a 

93 mixed R. ferrumequinum /M. emarginatus colony, in a parallel longitudinal study. 

94 Several novel viruses of both families were detected and we show that bats are also 

95 a host for BetaCoV type 1. 

96 MATERIALS 

97 In 2015 and 2016, fecal samples (n= 624) were collected from a mixed R. 

98 ferrumequinum /M. emarginatus nursing colony in Bech-Kleinmacher, using a 

99 longitudinal approach. Samples were collected (i) after resettling of the colony in the 

100 summer roost and before birth of the juveniles (June 2015, n=100; May 2016, n=99), 

101 (ii) during lactation (July 2015, n=126; June 2016, n=111) and (iii) before the colony 

102 returned to the winter roost (September 2015, n=100; September 2016, n=88). In 

103 201 6, in the framework of a cross-sectional study, fecal samples (n= 254) were 

104 collected from 6 of the 14 synanthropic M. emarginatus colonies known in 

105 Luxembourg (Table 1, Figure 1). Beginning of June 2016 and before birth of the 

106 juveniles, the population size of every known M. emarginatus maternity colony in 

107 Luxembourg (Table 1, Figure 1 A) was assessed by counting the bats emerging from 

108 the roost and/or the bats from a photograph taken in the roost, according to the 

109 Guidelines for Surveillance and Monitoring of European Bats (33). 

no The monitoring and sample collection was approved by the Ministry of Sustainable 

in Development and Infrastructure Luxembourg (ref.: 86503 CG/ne). 
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112 Fresh faeces were collected on a clean tarpaulin (left for 2-12 hours underneath the 

113 roost) and individually placed into 2 ml tubes using single-use spatulas. Samples 

114 were kept at +4°C during transport to the laboratory, where they were directly 

ns processed. Bat species was identified by visual inspection of the faeces and of the 

116 bat cluster hanging above the collection site. Species identification was confirmed for 

117 virus-positive samples by sequencing of mitochondrial DNA (see below). 

ns The study dataset is described in Table 1 and the primer sequences can be found in 

119 Table 2. 

120 METHODS 

121 Nucleic acid extraction. Entire bat droppings (approximately the size of a long grain 

122 of rice) were individually resuspended in 1 ml of prechilled virus transport medium 

123 (prepared according to the WHO protocol (34)) and homogenized using stainless 

124 steel beads (Qiagen, Venlo, The Netherlands) and a TissueLyser II (Qiagen). 

125 After centrifugation at 2200 g for 20 min, the supernatant was transferred to a new 2 

126 ml tube and stored at -80°C until further processing. Before nucleic acid extraction, 

127 each sample was centrifuged at 2200 g for 10 min and spiked with an extraction 

128 control (i.e. Human Adenovirus C5). Concurrent extraction of DNA and RNA was 

129 performed with the QIAamp Viral RNA Mini Kit (Qiagen) according to the 

130 manufacturer's protocol. To test for inhibition and confirm the successful extraction, 

131 each sample was tested using a real-time PCR specific for adenovirus (35). 

132 Virus detection. All samples were tested for CoV and PV by reverse transcription 

133 PCRs with degenerate primers in a nested format. The PCRs were performed in a 

134 final volume of 25 pi. In the first step of the nested PCR, the QIAGEN One-Step RT- 

135 PCR Kit (Qiagen) was used. The mastermix of the CoV PCR contained 2pl of RNA, 
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136 1 pM of each primer, 1 mM of MgCI 2 and 1 mM of each dNTP and the mastermix of 

137 the PV PCR contained 250nM of each primer, 1,5mM of MgCI 2 and 10OpM of each 

138 dNTP. In the second step of the nested PCRs, the mastermix of the CoV PCR 

139 contained 2.5 pi of 1:5 diluted PCR product, 700nM of each primer, 2mM of MgCI 2 

140 and 200pM of each dNTP, whereas the mastermix of the PV PCR contained 0.1 pi of 

141 undiluted PCR product, 600nM of each primer, 2mM of MgCI 2 and 200pM of each 

142 dNTP. The adenovirus detection PCR was similar to the CoV PCR, but 2.5 pi of DNA 

143 were used and 560nM of probe were added to the mix. In the second step of the 

144 nested PCRs, in the adenovirus detection PCR and in the bat species identification 

145 PCR, the Platinum® Taq DNA Polymerase Kit (Life Technologies Europe B.V., Gent, 

146 Belgium) was used. The CoV primers target the RNA-dependent RNA polymerase 

147 (RdRp) [modified from (36)], whereas the PV primers target the L-gene (37) of all 

148 known strains of the respective viral families. An avian infectious bronchitis virus (an 

149 avian CoV) and a measles virus (a human PV) served as positive controls in the CoV 

150 and PV PCRs. Details about the primers can be found in Table 2. 

151 Sequencing. PCR positive samples were identified by agarose gel-electrophoresis. 

152 Where multiple bands were present, amplicons of the appropriate size were excised 

153 from the gel and purified with the QIAquick gel extraction kit (Qiagen). PCR products 

154 giving a single band in the gel-electrophoresis were directly purified using the 

155 JetQuick™ extraction kit (Genomed, Lohne, Germany). Sequencing was performed 

156 using the BigDye terminator kit (Applied Biosystems, Foster City, CA) run on an ABI 

157 3130 sequencer (Applied Biosystems). Partial L-gene sequences of PV were 

158 obtained using the detection primers. Partial sequencing of CoV was attempted 

159 using specific primers targeting the conserved RdRp gene, as well as the spike 

160 glycoprotein gene. To reliably identify the bat species of all virus-positive samples, 
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161 partial cytochrome b sequences were obtained. The bat species identification PCR 

162 was performed using the Platinum® Taq DNA Polymerase Kit (Life Technologies 

163 Europe B.V.) in a final volume of 25 pi containing 5 pi of DNA, 700nM of each 

164 primer, 4mM of MgCI and 400pM of each dNTP. New primer sets were designed and 

165 evaluated with Geneious software (version 7.1.7; Biomatters Limited; Auckland, New 

166 Zealand [http://www.geneious.com]) and Primer3Plus 

167 (http://www.bioinformatics.nl/cgi-bin/primer3plus/primer3plus.cgi/). Details about the 

168 sequencing primers can be found in Table 2. 

169 All viral and mitochondrial sequences were submitted to GenBank (accession 

170 numbers KY502383 to KY502414, as well as KY707827 and MF048874 to 

171 MF048903). 

172 Sequence and phylogenetic analyses. Sequence assembling and processing was 

173 performed in Geneious v.7.1.9 (http://www.geneious.com/; (38)). A BLASTn search 

174 against the sequences in GenBank (https://www.ncbi.nlm.nih.gov/genbank/) was 

175 performed with the default parameters. Phylogenetic trees based on nucleotide 

176 sequences of the partial RdRp and spike genes for CoV, and of the partial L-gene for 

177 PV were constructed. In order to increase the phylogenetic resolution and because 

178 of the high genetic similarity of the obtained virus strains, only the longest sequences 

179 of high quality were selected among the novel sequences, and aligned with 

180 representative GenBank sequences using the ClustalW algorithm (39), as 

181 implemented in Geneious. Poorly aligned positions in the alignments were eliminated 

182 using Gblocks (40) as implemented in Seaview version 4 (41). Maximum likelihood 

183 (ML) and Bayesian inference of evolution were estimated in PhyML (42, 43) and 

184 BEAST (44, 45), respectively. The best substitution model identified by jModelTest 

185 (46), according to the Bayesian Information Criterion (BIC) and/or Akaike Information 
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186 Criterion (AIC) values was used. A bootstrap test including 1000 replicates was 

187 performed for each ML tree. For the Bayesian Markov Chain Monte Carlo (MCMC) 

188 approach, the parametric model “Constant Size” was used as prior and the analyses 

189 were performed with a lognormal relaxed clock. The MCMC run was at least 2X10 7 

190 steps long with sampling every 10 3 steps. Convergence was assessed on the basis 

191 of the effective sampling size using Tracer version 1.6. (47). The results of the 

192 Bayesian phylogenetic inference were summarized in a maximum clade credibility 

193 tree using the Tree Annotator program after a 10% burn-in. Tree topology was tested 

194 by posterior probability (pp) and only the pp values of well supported nodes (pp>0.7) 

195 are shown. As the topologies of the trees based on Bayesian and ML inference 

196 largely overlapped, only the maximum clade credibility trees are shown. However, for 

197 the nodes also supported by ML inference (bootstrap confidence levels above 0.7), 

198 the bootstrap support is shown in brackets. The scale bar of the trees indicates the 

199 average number of nucleotide substitutions per site (Figure 2-4). 

200 Statistical analyses. Statistical analyses were performed in R software (version 

201 3.1.0.; R Foundation for Statistical Computing, Vienna, Austria [https://www.r- 

202 project.org/]) (48). Logistic regression was performed to predict the binary outcome 

203 (i.e. presence or absence of detectable AlphaCoV shedding by M. emarginatus) 

204 based on the categorical predictor “season” with the levels “May/June”, “July” and 

205 “September” and using a logistic function. 

206 RESULTS 

207 To assess the prevalence and diversity of PV and CoV shedding among bats in 

208 Luxembourg, fecal samples from 7 colonies (Figure 1 A) were screened using 

209 degenerate primers in a nested format. The overall prevalence of PV was 1.1 % 

210 (10/878) and of CoV was 4.9% (43/878) and viruses were found in every colony 
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211 except for Colpach and Marienthal (Figure 1A, Table 1). No PV and CoV coinfections 

212 were detected. 

213 Bat PV were only detected in Ettelbruck and Bech-Kleinmacher (Figure 1 A), and 

214 shedding rates never exceeded 0.8-3.6% throughout the observation period. 

215 Because of the low prevalence rates, statistical analyses of the seasonal variation 

216 were not possible for PV. 9 out of the 10 detected PV strains were nearly identical to 

217 each other (represented by LUX15-A-033 and 351 in Figure 2). BLAST and 

218 phylogenetic analyses revealed that our PV strains were most closely related to 

219 those of insectivorous bats from China and South Africa with which they shared less 

220 than 80% nucleotide identity and less than 92% amino acid identity. Based on the 

221 phylogenetic analyses (Figure 2), all study sequences were grouped into a well- 

222 supported cluster, comprising also the unassigned murine J-virus (49), Beilong virus 

223 (50) and other Jeilongvirus-related PV (2, 3). 

224 From the CoV strains detected in this study, partial RdRp gene sequences were 

225 obtained. We show that 2 out of the 4 currently recognized CoV genera (i.e. 

226 AlphaCoV, BetaCoV, GammaCoV and DeltaCoV) circulate in Luxembourg and 36 

227 AlphaCoV and 7 BetaCoV were detected (Figure 3; Table 1). M. emarginatus from 

228 three different colonies (i.e. Ettelbruck, Lintgen and Bech-Kleinmacher; Figure 1 A, 

229 Table 1) shed nearly identical AlphaCoV{> 99 % nucleotide identity between partial 

230 RdRp gene sequences), most closely related to CoV circulating among insectivorous 

231 bats in China (Figure 3). In contrast to the PV shedding, AlphaCoV shedding was 

232 variable in Bech-Kleinmacher and highest rates were observed in July after 

233 parturition (OR: 2.8, p<0.01; Figure 1B). AlphaCoV strains from this study 

234 (represented by LUX15-A-48 in Figure 3) formed a distinct cluster and their RdRp 
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235 gene sequences shared less than 86% amino acid identity with other classified 

236 AlphaCo V. 

237 On the phylogenetic tree of the partial RdRp gene, the BetaCoV strains from this 

238 study clustered within 2 out of the 4 recognized lineages (A to D) of BetaCoV 

239 (https://talk.ictvonline.org): R. ferrumequinum shed SARS-related CoV (lineage B, 

240 represented by LUX16-A-24 in Figure 3) and M. emarginatus shed BetaCoV 1 

241 (lineage A, represented by LUX15-A-158 in Figure 3). The SARS-related CoV from 

242 Bech-Kleinmacher were identical to each other and BLAST analyses revealed 94% 

243 nucleotide identity between partial RdRp gene sequences of this study and SARS- 

244 related CoV circulating among rhinolophid bats in Europe (Figure 1A, Figure 3). 

245 Besides, we detected the first bat BetaCoV 1 strains (n= 5) in M. emarginatus from 3 

246 different colonies in 2015 (Bech-Kleinmacher) and 2016 (Bissen and Platen) (Figure 

247 1 A, Table 1). All strains from Luxembourg were highly similar to each other and to 

248 BetaCoV 1 identified in various mammalian species (>99% nucleotide identity 

249 between partial RdRp gene sequences) (Figures 3 and 4). 

250 Sequencing of the partial spike gene was attempted for all novel bat CoV strains, but 

251 was successful only for the BetaCoV 1 strains. As for the RdRp gene, all spike gene 

252 sequences were highly similar to each other and shared >98% nucleotide identity 

253 with the BetaCoV 1 strains from other mammalian species (Figure 4B). 

254 DISCUSSION 

255 Bats are natural reservoirs of numerous viruses with zoonotic potential. Of particular 

256 interest are CoV and PV that share several traits allowing their adaptation to new 

257 ecological niches and hosts: high mutation rates, poor RNA proofreading capability 

258 and genetic recombinations (51-53). In line with previous studies [e.g. (3, 4, 23)], we 
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259 found genetically diverse CoV and PV strains in bats that are known to forage in and 

260 around human settlements in Luxembourg (Figures 2 to 4, Table 1). Shedding rates 

261 may have been under-estimated due to RNA degradation, as well as due to low viral 

262 loads in faeces (54) and reduced sensitivity of degenerate primers. However, the 

263 sample collection and processing protocol was optimized to minimize the 

264 degradation of viral particle and of RNA, as well as inhibition. We acknowledge that 

265 the adenovirus control did not control for inhibition during the reverse transcription 

266 step. Although somewhat more susceptible to PCR inhibition and RNA degradation, 

267 fecal samples have been systematically used before to investigate virus 

268 epidemiology and evolution (55-58). Moreover, faeces are collected non-invasively 

269 and are thus the preferred material when studying viruses circulating among these 

270 endangered species (59, 60). 

271 Plowright et al proposed three scenarios to explain temporal variations in virus 

272 shedding in bats: (i) virus reactivation in persistently infected bats, (ii) seasonal 

273 epidemic cycles aligning with physiology of their life cycle, or (iii) transient epidemics 

274 due to waning immunity (19,61). Similar to a previous study (4), we observed no 

275 temporal variation of PV shedding possibly because of the low prevalence of PV 

276 shedding. In contrast and comparable to another study (55), significant increase in 

277 shedding of AlphaCoV was found in July possibly due to the periparturient stress 

278 (55, 62) (Figure IB). 

279 The lack of similar reference sequences complicated the genetic and phylogenetic 

280 characterization of the detected virus strains. Nevertheless, we identified novel PV 

281 and AlphaCoV strains that are related to bat viruses from distant regions of the world 

282 (Figures 2 and 3). Also, according to the PV species discrimination criterion 

283 published before (i.e. amino acid distance above 7-7.5% in the L-gene) (4), the study 
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284 sequences may represent putative novel PV strains, but this finding needs 

285 confirmation by whole genome sequencing. Based on the topologies of the 

286 phylogenetic trees and BLASTn analyses, the new CoV obtained in this study were 

287 found to be sufficiently divergent to represent a novel RGU based on amino acid 

288 sequence analysis of the partial RdRp gene (5, 23). We found no evidence of inter- 

289 species transmission, although a mixed-species colony was followed during 2 years 

290 (Figure 1). Together these findings confirm previous studies suggesting an 

291 association between AlphaCoV and host taxa rather than between geography and 

292 viral evolution, and thus a close virus-host evolution (23, 63-65). 

293 On the other hand, the detection of highly similar virus strains in different colonies 

294 (Figure 1 A, Table 1) is indicative of a social link between M. emarginatus colonies in 

295 Luxembourg. This is of particular interest with respect to ongoing efforts for the 

296 conservation of this species. Indeed, short foraging distances (24) and life-long roost 

297 fidelity complicate the preservation of M. emarginatus (66, 67). Since migratory 

298 distances of 35 to 126 km have been reported between summer and winter roosts 

299 (66, 67), and since all Luxembourgish colonies are within 45 km of each other 

300 (Figure 1 A), bats from different colonies may assemble during autumn swarming of 

301 the males (30, 68). Thus, male bats may play a particular role in virus transmission 

302 which warrants further investigation. Thus, a better understanding of the dynamics of 

303 bat-associated viruses may benefit these endangered species by indirectly providing 

304 information about foraging and mating behaviour. 

305 In contrast to AlphaCoV, host-switching is a major evolutionary mechanism of 

306 BetaCoV. For instance, SARS- and MERS-CoV circulated in bats, before crossing 

307 the species barrier to infect an intermediate host, which in turn infected humans (8, 

308 13, 69, 70). Bat SARS-CoVs even use the same receptor for cell entry than their 
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309 human counterpart and they were detected in rhinolophid bats (8) that host also 

3 10 genetically diverse SARS-related CoV (23, 71-74). Also in our study, Ft. 

311 ferrumequinum from Bech-Kleinmacher shed SARS-related CoV strains (Figure 3). 

312 Although it is unlikely that these CoV represent a direct threat for humans, the 

313 potential risk of adaptation to the human host should not be ignored (75-77). 

314 BetaCoV 1 species is another exception to the typical host-specificity of CoV. The 

315 species comprises highly similar viruses of distantly related mammals (6, 78-81) and 

316 so far, only a single, short BetaCoV 1 sequence was obtained from a bat (10). Most 

317 interestingly, we show here that M. emarginatus from different roosts shed BetaCoV 

318 1 strains (Figure 1 A) that are highly similar and closely related to BetaCoV 1 strains 

319 detected in various other amimal species (Figure 4). Most recent common ancestor 

320 analyses of BetaCoV 1 suggested that the group appeared only recently and has low 

321 host-specificity (82-84). For example, BetaCoV 1 strains detected in exotic ruminants 

322 such as giraffes or antelopes are thought to represent spillover viruses of bovine 

323 CoV that underwent adaptive mutations (78, 80). Moreover, a possible animal origin 

324 of the human FICoV OC43 was revealed by molecular clock analysis of the spike 

325 gene (83, 84) that provides indication about host range and tissue tropism. The 

326 permissiveness of human cells to certain BetaCoV 1 strains further underlines the 

327 potential of BetaCoV 1 to be transmitted across species (80, 82). Also here, all spike 

328 gene sequences were highly similar to BetaCoV 1 from other mammalian species 

329 reflecting the genetic stability typical of the lineage (78, 85, 86). To further investigate 

330 the role of bats as reservoir of BetaCoV 1, studies focusing on the host range of this 

331 CoV species are warranted. 

332 In conclusion, we showed that bats in Luxembourg, Western Europe, are hosts of 

333 novel virus strains that may be able to overcome the species barrier. BetaCoV 1 
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334 strains with spike and RdRp genes genetically highly similar to mammalian strains 

335 were detected in synanthropic bats. In addition, we identified SARS-related CoV that 

336 may infect humans after a viral adaptation process (75-77). As shown before for bat 

337 lyssaviruses (32), our study highlights a certain risk for zoonotic transmission of bat 

338 viruses in particular since foraging and roosting sites of most indigenous bat species 

339 overlap with human and animal habitats. To mitigate this risk, it is important to 

340 monitor viruses circulating in synanthropic bats and putative intermediate hosts, and 

341 to identify factors that affect bat populations. 
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613 

614 

TABLES AND FIGURE LEGENDS 

Table 1 Characteristics of the different colonies and of the dataset, as well as detection rates of coronaviruses (CoV) and 

paramyxoviruses (PV) 

Location of Bat species Population Sample Sample Viral nucleic acids detected in bat faeces 

colony (UTM size on the collection size (% allCoV* ana ** AlphaCoV* SARS- BetaCoV 1* bat PV* (%) 

coord.) 10/06/2016 date subset of (%) (%) related (%) 

population) CoV** (%) 



Lintgen (32U 

Myotis 

65 

10.06.2016 

44 (67.7) 

8 (18.2) 

8 (18.2) 

0(0) 

0(0) 

0(0) 



293/5511) 

emarginatus 









o 


Ettelbruck 

Myotis 

220 

10.06.2016 

44 (20) 

1 (2.3) 

1 (2.3) 

0(0) 

0(0) 

1 (2.3) 

c 


(32U 

emarginatus 









£ 


291/5525) 












Marienthal 

Myotis 

45 

03.06.2016 

45(100) 

0(0) 

0(0) 

0(0) 

0(0) 

0(0) 

P CD 


(32U 

emarginatus 









> “5 


288/5510) 










c .2 


Bissen (32U 

Myotis 

35 

03.06.2016 

33 (94.3) 

1 (3) 

0(0) 

0(0) 

1 (3) 

0(0) 



288/5519) 

emarginatus 









c b 


Colpach 

Myotis 

70 

10.06.2016 

44 (62.9) 

0(0) 

0(0) 

0(0) 

0(0) 

0(0) 

o 


(31U 

emarginatus 









~u ^ 


703/5515) 












Platen (31U 

Myotis 

60 

16.06.2016 

44 (73.3) 

2 (4.6) 

0(0) 

0(0) 

2 (4.6) 

0(0) 

Q_ 


711/5519) 

emarginatus 











Bech- 

Myotis 

942 

09.06.2015 

100(10.6) 

1 (1) 

0(0) 

0(0) 

1 (D 

1 (1) 



Kleinmacher 

emarginatus 


14.07.2015 

126(13.4) 

15(11.9) 

15 (11.9) 

0(0) 

0(0) 

1 (0.8) 



(32U 

and 


04.09.2015 

100(10.6) 

7(7) 

6(6) 

0(0) 

1 (D 

2(2) 



308/5489) 

Rhinolophus 


17.05.2016 

99(10.5) 

3(3) 

1 (1) 

2(2) 

0(0) 

1 (1) 




ferrumequinum 

11.07.2016 

111 (11.8) 

5 (4.5) 

5 (4.5) 

0(0) 

0(0) 

4 (3.6) 






09.09.2016 

88 (9.3) 

0(0) 

0(0) 

0(0) 

0(0) 

0(0) 



Total 


1437 


878 (61.1) 

43 (4.9) 

36 (4.1) 

2 (0.2) 

5 (0.6) 

10(1.1) 


615 

* detected in Myotis emarginatus , ** i 

detected in Rhinolophus ferrumequinus, *** UTM coord.: Universal Transverse Mercator 


616 

coordinates 
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618 Table 2 Primers used for detection and sequencing of Coronaviruses and Paramyxoviruses 


Virus 

Tarqet protein 

Primer sense 

Primer sequence (5'-3') 

Amplicon size (bp) 

Reference 

Detection 






Paramyxoviridae 

L protein, subunit of 

Forward 

G AAGG IT ATT GT CAIAARNTNT GGAC 

660 

(37) 


RNA-dependent RNA 

Reverse 

GCTGAAGTTACIGGITCICCDATRTTNC 




polymerase 

Forward semi-nested 

GTTGCTT CAAT GGTT CARGGNGAYAA 

580 




PCR 




Coronaviridae 

Replicase polyprotein 

Forward 

GGKTGGGAYTAYCCKAARTG 

602 

(36) 


lab 

Reverse 

TG YTGT SW RCARAAYTCRT G 





Forward nested PCR 

GGTTGGG ACT AT CCT AAGT GT G A 

555 




Reverse nested PCR 

CCAACAYTTNGART CWGCCAT 


this study 

Human 

Hexon gene 

Forward 

GCCACSGTGGGGTTYCT AAACTT 

130 

this study 

adenovirus 








Reverse 

GCCSC AGT GG KCDT AC AT GCACAT C 


this study 



Probe 

FAM-TGCACCAG ACCCGGGCT CAGGT ACT CCG A-T AM R A 


(35) 

Sequencing 






Alphacoronavirus 

Replicase polyprotein 

Forward 

TGTGAAGGCCTTACAGCGTC 

670 

this study 


lab 

Reverse 

AGAGCCACAWACAACACACA 




Replicase polyprotein 

Forward 

T GAT GCAGCT GTYARAGACTTC 

690 

this study 


lab 

Reverse 

CCAGAAGTCGTACCACCAGG 



Betacoronavirus 

Replicase polyprotein 

Forward 

AGACATCGTCCCCATCCATC 

729 

this study 

1 

lab 

Reverse 

AGCT ACACGTGGT GTT CCT G 




Replicase polyprotein 

Forward 

CAT AT CAT CCCAGCCGCCAT 

584 

this study 


lab 

Reverse 

TG CT GTTTT AGT GTT G CGG C 




Replicase polyprotein 

Forward 

CCGCTT GTT AT AGCCGCAAC 

613 

this study 


lab 

Reverse 

AGCGCT ACT GAGTTTGCAGA 




Spike gene 

Forward 

GTGAGCACTGTTCGGGTCTT 

432 

this study 



Reverse 

AGCAAT GCTGGTT CGGAAG A 




Spike gene 

Forward 

AT GGC ATT GGG AT ACAG 

492 

(80) 



Reverse 

T AATGGAGAGGGCACCGACTT 




Spike gene 

Forward 

GGGTT ACACCT CT CACTT CT 

767 

(80) 



Reverse 

GCAGGACAAGTGCCTATACC 



SARS-related 

Replicase polyprotein 

Forward 

AGTT GAGGTGGTCGACAAGT 

650 

this study 

Coronavirus 

lab 

Reverse 

GCAGT GGT AGCAT CT CCT GA 



Bat species 

Cytochrome b 

Forward 

ATGACCAACATTCGMAARTCYCAC 

390 

this study 

identification 








Reverse 

T GAT GACGGTTGCTCCTCA 
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619 Figure 1 Municipalities with known Myotis emarginatus colonies in 

620 Luxembourg and circulation of coronavirus and paramyxovirus strains (A); 

621 Seasonality of Alphacoronavirus shedding in Bech-Kleinmacher (B). 

622 (A) Sampled municipalities are in bold. The blue quadrant with the mixed colony of 

623 Bech-Kleinmacher as centre, has a radius of 45 km and embraces all investigated 

624 colonies. The base map is from the Land Registry Office of the Grand Duchy of 

625 Luxembourg. (B) Error bars represent the 95% confidence interval and * 

626 corresponds to p<0.05. 

627 Figure 2 Phylogenetic analysis of the partial L gene of paramyxoviridae. 

628 Bayesian analyses of a 410 nt long alignment comprising unique, partial L gene 

629 sequences of 34 paramyxovirus (PV) strains that represent all PV species 

630 recognized by the International Committee on Taxonomy of Viruses (ICTV), as well 

631 as unassigned, but putative novel PV species. Three out of the 10 PV from this study 

632 were added to the dataset to represent the genetic diversity of PV circulating in 

633 Myotis emarginatus in Luxembourg. Four Pneumoviridae strains served as outgroup 

634 for the phylogenetic analyses. The study sequences are in red and strains hosted by 

635 bats are highlighted in bold to show the high genetic diversity of bat PV. Only the pp 

636 values of well supported nodes (pp>0.7) are shown and if the nodes were also 

637 supported by ML inference (bootstrap confidence levels above 0.7), the bootstrap 

638 support is shown in brackets. For each cluster, the PV species, as well as the virus 

639 family assignment are shown. The sequences were named, if the information was 

640 available, according to the following nomenclature: abbreviated virus name/host 

641 species/three-letter code of the country of origin/Genbank accession number. The 

642 following virus name abbreviations were used: PMPV-Pneumonia virus of mice; 
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643 HRSV-Human respiratory syncytial virus; HMPV-Human metapneumovirus; AMPV- 

644 Avian metapneumovirus type; HPIV-Human parainfluenza virus; MapV-Mapuera 

645 virus; LPMV- Porcine rubulavirus; SV-Simian virus; MuV-Mumps virus; TuV-Tuhoko 

646 virus: TV-Tupaia paramyxovirus; MenV-Menangle virus; APMV-Avian 

647 paramaxyovirus; MeV-Measles virus; PPRV-Peste des petits ruminants virus; DoV- 

648 Dolphin morbillivirus; PDV-Phocine distemper virus; CDV-Canine distemper virus; 

649 NaV-Nariva virus; MoV Mossman virus; BtV-Bat paramyxovirus; JV-J-virus; BV- 

650 Beilong virus; NiV-Nipah virus; HeV-Hendra virus; BPIV-Bovine parainfluenza virus; 

651 PPIV-Swine parainfluenza virus; HPIV-Human parainfluenza virus; SeV-Sendai 

652 virus; AsaPV- Atlantic salmon paramyxovirus: AnV-Anaconda paramaxyovirus; 

653 FDLV- Fer-de-Lance paramyxovirus. 

654 Figure 3 Phylogenetic analysis of the partial RNA-dependent RNA polymerase 

655 gene of all Coronavirus genera 

656 Bayesian analyses of a 853 nt long alignment comprising unique, partial RNA- 

657 dependent RNA polymerase gene sequences of 50 coronavirus (CoV) strains that 

658 represent all CoV species recognized by the International Committee on Taxonomy 

659 of Viruses (ICTV), as well as unassigned, but putative novel CoV species. Three out 

660 of the 43 CoV from this study were added to the dataset to represent the genetic 

661 diversity of CoV circulating in Myotis emarginatus and Rhinolophus ferrumequinum 

662 in Luxembourg. The DeltaCoV strains served as outgroup for the phylogenetic 

663 analyses. Only the pp values of well supported nodes (pp>0.7) are shown and if the 

664 nodes were also supported by ML inference (bootstrap confidence levels above 0.7), 

665 the bootstrap support is shown in brackets. 
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666 For each strain, the CoV genus assignment is shown. Assignment to recognized 

667 ICTV species is only shown for the study sequences that are highlighted in red. The 

668 country of origin of each strain is stressed in bold to show the vast geographic 

669 spread of CoV. The sequences were named, if the information was available, 

670 according to the following nomenclature: abbreviated virus name/virus strain/host 

671 species/three-letter code of the country of origin/year of sampling/specimen number^ 

672 Genbank accession number. The following virus name abbreviations were used: 

673 BtCoV-Bat coronavirus; PEDV-Porcine epidemic diarrhea virus; HCoV-Human 

674 coronavirus; TGEV- Transmissible gastroenteritis virus; PRCV-Porcine respiratory 

675 coronavirus; FCoV-Feline coronavirus; MiCoV-Mink coronavirus; CiCoV-Civet 

676 Severe acute respiratory syndrome CoV; RatCoV-Rat coronavirus; MFIV-Murine 

677 hepatitis virus; PFIEV-Porcine hemagglutinating encephalomyelitis virus; BCoV- 

678 Bovine coronavirus; CaCoV-Canine respiratory coronavirus; ECoV-Equine 

679 coronavirus; DrCoV- Dromedary camel coronavirus; WtDCoV-White-tailed deer 

680 coronavirus; GCoV-Giraffe coronavirus; AnCoV-Sable antelope coronavirus; 

681 WBkCoV-Waterbuck coronavirus; SdCoV-Sambar deer coronavirus; RabCoV-Rabbit 

682 coronavirus; HeCoV-Hedgehog coronavirus; IBV-Infectious bronchitis virus; CMCoV- 

683 Common moorhen coronavirus; WECoV-Wigeon coronavirus; BuCoV-Bulbul 

684 coronavirus; ThCoV-Thrush coronavirus; MuCoV-Munia coronavirus; PCoV-Porcine 

685 coronavirus; WiCoV-White-eye coronavirus; NFICoV-Night heron coronavirus. 

686 Figure 4 Phylogenetic analysis of the partial RNA-dependent RNA polymerase 

687 gene (A) and partial spike glykoprotein gene (B) of Betacoronaviruses 

688 Bayesian analyses of a 1771 nt long alignment comprising unique, partial RNA- 

689 dependent RNA polymerase gene sequences of 38 coronavirus (CoV) strains (A) 

690 and of a 911 nt long alignment comprising unique, partial spike glykoprotein gene 
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691 sequences of 34 CoV strains that represent all BetaCoV species recognized by the 

692 International Committee on Taxonomy of Viruses (ICTV). In addition, 1 out of the 5 

693 highly similar BetaCoV 1 strains from this study was added to the dataset to show 

694 the genetic relationship of BetaCoV 1 circulating in Myotis emarginatus in 

695 Luxembourg to BetaCoV 1 of other host species. Only the pp values of well 

696 supported nodes (pp>0.7) are shown and if the nodes were also supported by ML 

697 inference (bootstrap confidence levels above 0.7), the bootstrap support is shown in 

698 brackets. Assignment to recognized ICTV species is shown for each strain. The 

699 study sequence is highlighted in red and strains that were detected in bats are 

700 stressed in bold to show that most BetaCoV species comprise CoV strains that were 

701 initially detected in bats. The sequences were named, if the information was 

702 available, according to the following nomenclature: abbreviated virus name/virus 

703 strain/host species/three-letter code of the country of origin/year of 

704 sampling/specimen number_Genbank accession number. The following virus name 

705 abbreviations were used: BtCoV-Bat coronavirus; HCoV-Human coronavirus; CiCoV- 

706 Civet Severe acute respiratory syndrome CoV; RatCoV-Rat coronavirus; MHV- 

707 Murine hepatitis virus; PHEV-Porcine hemagglutinating encephalomyelitis virus; 

708 BCoV-Bovine coronavirus; CaCoV-Canine respiratory coronavirus; ECoV-Equine 

709 coronavirus; DrCoV- Dromedary camel coronavirus; WtDCoV-White-tailed deer 

7 10 coronavirus; GCoV-Giraffe coronavirus; AnCoV-Sable antelope coronavirus; 

711 WBkCoV-Waterbuck coronavirus; SdCoV-Sambar deer coronavirus; RabCoV-Rabbit 

712 coronavirus; HeCoV-Hedgehog coronavirus. 
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PMPV/mouse/USA_AY743910 
HRSV-A2/human_HRU39662 
HMPV/human/NLD_AY525843 
AMPV/bird/USA_AY394492 
HPIV-4b/human/CAN_EU627591 
MapV/fruit bat/BRA_EF095490 
LPMV/pig/MEX_BK005918 
SV-41_NC006428 
MuV/human/RUS_AY669145 
TuV-1/fruit bat/CHN_GU128080 
TiV/fruit bat/MYS_AF298895 
MenV/fruit bat/AUS_AF326114 
TuV-3/fruit bat/CHN_GU 128082 
TuV-2/fruit bat/CHN_GU 128081 ? 
APMV-9/duck/USA_EU910942 A 

APMV-1/bird/AUS_AY935496 ▼ 

MeV/human/JPN_AB254456 
PPRV/sheep/TUR_AJ849636 
DoV/dolphin_AJ608288 
PDV/seal_Y09630 
CDV/raccoon_AY443350 
NaV/wild rodent/TTO_FJ362497 
MoV/wild rat/AUS_AY286409 
BtV/bat/DEU_JN086950 
LUX15_A_033/bat/LUX_KY502406 
LUX15_A_351 /bat/LUX_KY502409 
LUX15_A_130/bat/LUX_KY502407 
BtV/bat/ZAF_KC578570 
JV/Mus musculus/AUS_NC007454 
BV/bat_DQ1 00461 
NiV/pig/MYS_AJ564623 
HeV/AUS_AF017149 
BtV/bat/ZMB_AB853096 
TV/rodent/THA_AF079780 
BPIV-3/cattle/USA_AF1 78654 
PPIV-1/pig/USA_EU439428 
HPIV-1/human/USA_NC003461 
SeV/rodent/JPN_AB039658 
AsaPV/salmon/NOR_EU 156171 
AnV/reptile/HKG_KJ956405 
FDLV/reptile/_AY 141760 


PV SPECIES 

^OrthopneumoV 

0MetapneumoV 


lunassigned PV 


■JeilongV-related 


^HenipaV 
0unassigned PV 

^RespiroV 

0AquaPV 

0FerlaV 


VIRUS FAMILY 

Pneumoviridae 


Paramyxoviridae 
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CoV SPECIES 


CoV GENUS 


0.7 


1 


jm 



BtCoV/HKU7-1/bat_DQ249226 
BtCoV/HKU8/Miniopterus bat/HKG_NC010438 
BtCoV/AFCD307/bat/HKG_EU420137 
BtCoV/AFCD62/bat/HKG_NC010437 
PEDV/CV777/pig_NC003436 

BtCoV/CDPHEI5/Myotis lucifugus/USA/2006_NC022103 
BtCoV/512/Scotophilus bat/CHN/2005_NC009657 
BtCoV/LUXI 5_A_48/Myotis emarginatus/LUX/2015_KY502384 
BtCoV/HKUlO/Rousettus bat/CHN/2005_NC018871 
HCoV/NL63/human_NC005831 
HCoV/229E/human_NC002645 
BtCoV/HKU2-33/Rhinolophus bat/HKG/2006_EF203067 
TGEV/SHXB/Sus scrofa/CHN/2013_KP202848 
PRCV/OH7269/pig/USA/2014_KR270796 
FCoV/DF-2 R3i/Felis catus_JQ408980 
MiCoV/WD 1127/Mustela vison/USA/1998_NC023760 
BtCoV/LUXI 6_A_24/Rhinolophus ferrumequinum/LUX/2016_KY502395, 
BtCoV/Rfl/Rhinolophus ferrumequinum/2004_DQ412042 

- - - 1 AY463060 ' 


HCoV/ShanghaiQXC2/human/CHN_, 
CiCoV/SZI 6/civet/HKG AY304488 


0 unassigned CoV 


SARS-related CoV 


0.94 


0.99 

m 




RatCoV/Parker/rat_NC012936 
MHV/Murine hepatitis virus/mouse_AF029248 
HCoV/HKU1/human/HKG_DQ415901 
PHEV/pig/BEL_DQ011855 
BCoV/Mebus/cow_U00735 
CaCoV/K37/dog/KOR/2008_JX860640 
oWCoV/4408/human/DEU/1988_FJ415324 
ECoV/NC99/horse_E F446615 

BtCoV/LUXI 5_A_158/Myotis emarginatus/LUX/2015_KY502388 
DrCoV/HKU23-362F/dromedary/ARE/2013_KF906250 
WtDCoV/OH-WD470/white-tailed deer/USA/1994_FJ425187 
GCoV/OH3-TC/giraffe/USA/2003_EF424622 
AnCoV/OH 1 /antelope/USA/2003_EF424621 
WBkCoV/OH-WD358/waterbuck/USA/1994_FJ425186 
SdCoV/OH-WD388-TC/sambar deer/USA/1994_FJ425190 
HCoV/OC43/human/USA AY585228 


RabCoV/HKU14-1/Oryctolagus cuniculus/CHN/2006_JN874559 

-RatCoV/HKU24/Rattus norvegicus/CHN/2012_KM349743 

■ BtCoV/HKU9-1/bat/CHN_EF065513 




BetaCoV 1 





HeCoVA/MC/Erinaceus europaeus/DEU/2012_KC545386 
BtCoV/H KU5-5/bat/C H N_EF065512 
HCoV/MERS 2874/human/OMN/2013_KT156561 
BtCoV/SC2013/Vespertilio superans/CHN/2013_KJ473821 
BtCoV/HKU4-4/bat/CHN_EF065508 
IBV/IBV-p65/Beaudette_DQ001339 
CMCoV/HKU21/common moorhen/HKG/2007_JQ065049 
WECoV/HKU16-6847/white-eye/HKG/2007_JQ065044 
BuCoV/HKU 11 -934/red-whiskered bulbul/HKG/2007_FJ376619 
ThCoV/HKUl 2-600/grey-backed thrush/HKG/2007_NC011549 
MuCoV/HKUl 3-3514/white-rumped munia/HKG/2007_NC011550 
PCoV/HKUl 5-155/pig/H KG/2010_JQ065043 
WiCoV/HKU20-9243/wigeon/HKG/2008_JQ065048 
NHCoV/HKU 19-6918/night-heron/HKG/2007_JQ065047 


AlphaCoV 


BetaCoV 

| Gamma CoV 
I DeltaCoV 


0.2 
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A 



I 


CoV Species 


Severe acute respiratory 
syndrome (SARS)-related CoV 


unassigned 

unassigned 

unassigned 


♦ 

♦ 


Rousettus bat CoV HKU9 
HumanCoVHKUl 

Murine CoV 1 


Betacoronavirus 1 




♦ 


Hedgehog CoV 1 
Pipistrellus bat CoV HKU5 
Tylonycteris bat CoV HKU4 

Middle East respiratory 
syndrome(MERS)-related CoV 


B 





:CoV/Rf1/Rhinolophusferrumequinum/2004_DQ412042 


■ r 

I iwifci 


.BtCoV/HKU3-1/bat/CHN_DQ022305 
.BtCoV/BM48-31/Rhinolophus blasii/BGR/2008_NC_014470 
,CoV/ShanghaiQXC2/human/CHN_AY463060 
.CiCoV/SZI 6/ci vet/H KG_AY304488 
,BtCoV/WIV1/Rhinolophus sinicus/CHN/2012_KF367457 
“ ':oV/Zhejiang2013/Hipposideros pratti/CHN/2013_KF636752 



BtCoV/HKU5-1/bat/CHN_EF065509 
:CoV/SC2013/Vespertilio superans/CH N/2013_K J473821 
HCoV/MERS 2874/human/OMN/2013_KT156561 
BtCoV/HKU4-1/bat/CHN_EF065505 
BtCoV/PML-PHEI/Neoromicia capensis/ZAF/2011_KC869678 
,HeCoVA/MC/Erinaceus europaeus/DEU/2012_KC545386 
, RatCoV/Pa rker/rat_N C012936 
,HCoV/HKU1/human/HKG_DQ415901 
.MHV/HKU24/Rattus norvegicus/CHN/2012_KM349743 
,ECoV/NC99/horse_EF446615 
,PHEV/PHEV/pig/BEL_DQ011855 
■ iCoV/OH1/antelope/USA/2003_EF424621 
lGCoV/OH3-TC/giraffe/USA/2003_EF424622 
BCoV/AH 187/cow/USA/2000_FJ938065 
,WBkCoV/OH-WD358/waterbuck/USA/1994_FJ425186 
lsdCoV/OH-WD388-TC/sambar deer/USA/1994_FJ425190 
,WtDCoV/OH-WD470/white-tailed deer/USA/1994_FJ425187 
i BCoV/Mebus/cow_U00735 
,HCoV/4408/human/DEU/1988_FJ415324 
,BtCoV/LUX15_A_158/Myotis emarginatus/LUX/2015_KY502388 
,HCoV/OC43/human/USA_AY585228 
iCaCoV/K37/dog/KOR/2008_JX860640 
,RDrCoV/HKU23-362F/dromedary camel/ARE/2013_KF906250 
RabCoV/HKU14-1/Oryctolagus cuniculus/CHN/2006_JN874559 
BtCoV/BtKY24/Eidolon helvum/KEN/2006_HQ728482 
BtCoV/HKU9-4/bat/CHN_EF065516 
BtCoV/HKU9-3/bat/CH N_EF065515 
BtCoV/HKU9-2/bat/CH N EF065514 


CoV Species 

^ SARS-related CoV 
0 unassigned 
^ SARS-related CoV 

0 unassigned 
• Pipistrellus bat CoV HKUS 

^ MERS-related CoV 
0 Tylonycteris bat CoV HKU4 
0 MERS-related CoV 
0 Hedgehog CoV 1 
0 Murine CoV 1 
0 Human CoVHKUl 
0 Murine CoV 1 


Betacoronavirus 1 


0 unassigned 
A Rousettus bat CoV HKU9 
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